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1. Introduction 

In this paper I present a brief survey of the active area of Special Functions associated with 
Root Systems. The article is intended for a general mathematical audience. It will not 
suppose prerequisites on either special functions or root systems. It will also skip many 
technical details. Some early work in this area (root systems BC2 and A2) was done by 
the author [17] and Sprinkhuizen [32]. During the last ten years important break-throughs 
were made by Heckman and Opdam [11], [27], [28], [12], Macdonald [23], [24], Dunkl [8] 
and Cherednik [3], [4], [5]. 

A lot of the motivation for the subject of this paper comes from analysis on semisimple 
Lie groups. Spherical functions on Riemannian symmetric spaces of the compact or non- 
compact type can be written as special functions depending on parameters which assume 
only special discrete values. In the one- variable cases these special functions are classical, 
also for parameter values without group theoretic interpretation, but in the more- variable 
cases they were new. In the case of group theoretic interpretation, many properties of 
these special functions, as well as associated harmonic analysis, immediately follow by 
group theoretic arguments. The case of more general parameter values yields the special 
functions associated with root systems. Properties derived in the group case can still be 
formulated in the general case, but now as conjectures rather than theorems. This paper 
describes some of the progress which has been made in proving these conjectures. For 
convenience, I will restrict to the polynomial (compact) case, with Bessel functions as a 
sole exception. (For an introduction to the non-polynomial case see [29], [13].) Neither will 
I discuss the recent work on commuting operators with elliptic functions as coefficients. 
An important aspect of the whole theory, which will not be discussed very much in this 
paper, is the connection with completely integrable systems, for instance the generalized 
Calogero-Moser system. 

Special functions associated with root systems have also been developed in the g-case, 
where g is a deformation parameter giving back the earlier cases when q = 1. Motivation 
and development of the theory in the g-case has been quite different from the q = 1 case. 
Except for the case of Hall polynomials [22], theory was developed [23], [24], [18] without 
interpretation in group theory. But afterwards quantum groups looked very promising as 
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a natural setting for these polynomials. This had turned out to be true in the one- variable 
case [19], and very recently some interpretations of more- variable cases on quantum groups 
were found [25], [10]. 

In any case, a quantum group interpretation for generic values of the parameters can- 
not be expected. But, by Cherednik's work [3], [4], [5] we know already another algebraic 
setting for special functions associated with root systems: afRne and graded Hecke algebras 
[20]. As shown by work of Opdam [29], this new algebraic context also allows harmonic 
analysis. 



2. The one- variable case 

In this section I will introduce three classical families of special functions, each depending 
on a real parameter A; > 0, and such that the cases k = and k = 1 are elementary. The 
three families are connected with each other by limit transitions. Later, for each of the 
families I will discuss generalizations which are associated with root systems. 

2.1. Bessel functions. Consider Bessel functions in a non-standard notation: 

Here we use the notation for shifted factorial: 

{a)j ■=a{a + l)...{a + j -1) (j = l,2, ...); (a)o := 1. 

The function J^k is related to the Bessel function Jq, in standard notation [9, Ch. VII] by 



Jk{x) = — T —Jk-i{x). 

X 2 2 

Note that 

Jk{x) = Jk{-x), Jfc(O) = 1. (2.2) 

The cases k = and k = 1 yield elementary functions: 

Jo{x) = cosx, Ji{x) = ^^^ . (2.3) 

X 

The function x \-^ Jk{\x) (A G M) is eigenfunction of a differential operator: 

d? 2k d ^ _ ,, . ,2 



It is the unique C°° solution of this differential equation under conditions (2.2). 



2.2. Ultraspherical polynomials. Consider ultraspherical or Gegenbauer polynomials 
[9, §10.9], i.e. polynomials C^ of degree n on R such that 

/ C^(cosa;)C^(cosa;) (sina;)^''(ia: = {n,m E Z+, n ^ m). 
Jo 

Then the C„ are determined up to a constant factor (in general, we will not use the 
standard normalization for Gegenbauer polynomials). For A; = 0, 1 we have: 

Q 1 sin((n + l)a;) 

C^fcosa;) = const, cosfnx), C„{cosx) = const. . (2.4) 

sma; 

The function x h^ C^(cosa;) is eigenfunction of a differential operator: 



dx^ 



2k cot a;— I C^(cosa;) = -n{n + 2k)C^{cosx). 



For (njv) being a sequence of positive integers such that un/N -^ A for some A > as 
iV ^ oo, we have the limit result 

C!^^(cos(x/N)) , ^ 

2.3. Q-Ultraspherical polynomials. Let < q < 1 and define for any complex a: 

CXD 

(a;g)oo := JJ(1 - aq^). 

j=0 

The infinite product converges because of the condition on q. We will consider q-ultras- 
plierical polynomials [1] in a non-standard notation. These are polynomials C^'*^ of degree 
ri on R such that 



^C^<^(cosa.)C^'^(cosa.) 

^0 



,2ix. ^ 2 



(e^-;g). 



dx = (n, m G Z_|_, n ^ m). 



Then the C^'*^ are determined up to a constant factor. If we put 

P„(e^^):=C^^(cosa:) 
then Pn is eigenfunction of a q-difference operator: 

Note that the P„ on the left hand side have arguments off the unit circle, while orthog- 
onality is on the unit circle. The cases k = and k = 1 are elementary as in (2.4) (not 
depending on q): 

Q ^1 a/ \ sin((n + l)a:;) 

C„'^(cosa;) = const. cos(na;), C„'^(cosar) = const. -. . 

sma? 

With suitable normalization there is the limit relation 

limC^''^(cosa;) = C^(cosa;). 

The g-ultraspherical polynomials form a subclass of the Askey-Wilson polynomials [2]: a 
family of orthogonal polynomials depending, apart from q, on four non-trivial parameters. 



2.4. Dunkl operators in one variable. We will now generalize the elementary for- 
mulas 

^t\x ^ j^^^Xx) + i\x J^{\x) and — e^^^ = i\ e^^^ (2.5) 

ax 

(the first formula follows by (2.3)). Dunkl [8] generalized the operator d/dx to a mixture 
of a differential and a reflection operator: 

{D^^^f){x) := fix) + k -^^^^ ~ ■^^""'^ . (2.6) 

X 

Note that this Dunkl operator sends smooth functions to smooth functions. Let us define 
a generalized exponential function in terms of Bessel functions (2.1) by 

£k{Xx) := Jk{Xx) + -f^ Jk+ii\x). (2.7) 

Then it follows immediately from well-known differential recurrence formulas for Bessel 
functions that 

D'-''^£k{Xx) = iX£k{Xx). (2.8) 

Formulas (2.7) and (2.8) generalize the formulas in (2.5). The function x h^ £k{Xx) is the 
unique C°° function which equals 1 in and which is eigenfunction with eigenvalue iX of 

For (Z)^''))^ we compute 

{D^'^? fix) = fix) + ^ fix) - k ^^"^"/^""^ 

X X'^ 

Thus, on even functions / the square of the Dunkl operator acts as the differential operator 
{d/dxY + 2kx~^ d/dx. In particular, its action on 

Jk{Xx) = ^{£k{Xx) + £k{-Xx)) 

yields 

{D^''^fJk{Xx) = -X^Jk{Xx). 



3. Preliminaries about root systems 

3.1. Definition of root system. Let V be a d-dimensional real vector space with inner 
product ( . , . ). For a. G ^\{0} let Sq, denote the orthogonal reflection with respect to the 
hyperplane orthogonal to a. (cf. Fig. 1): 

.,(/3):=/5-^^a (/3 G F). 



a 



,a) 



A root system in V (see [15]) is a finite subset R of y\{0} which spans V and which 
satisfies for all a, (3 E R the two properties that 

s^{(3) G it! and / ' / G Z. 

(a, a) 

Clearly, if a G i? then —a = Sa{a) G -R. For convenience, we will restrict ourselves to the 
case of a reduced root system, i.e., a root system R such that, if a,/5 G i? and a = c/S 
for some c G M, then c = ±1. The so-called irreducible root systems can be classified 
as four infinite families An, Bn, Cn, Dn of classical root systems and five exceptional root 
systems G2,F4,Eq, E-j, Eg. Here the subscript denotes the rank of the root system, i.e. the 
dimension of V. There is one infinite family of non-reduced irreducible root systems: of 
type BCn- 




^ 
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Fig. 1. Refiection 5^ 



Fig. 2. Root system Ai 



Fig. 3. Root system C2 



An example for d = 1 is the set R := {±2} C K (root system of type Ai, cf. Fig. 2). An 
example for d = 2 is the set R = R+U {-R+), where R+ := {(1, -1), (2, 0), (1, 1), (0, 2)} C 
M? (root system of type C2, cf. Fig. 3). In general, when we have a root system R in 
V then we can write it as a disjoint union R = i?_|_ U (— i?_|_), where i?_|_ and —R+ are 
separated from each other by a hyperplane in V through the origin. The choice for i?_|_ is 
not unique. The elements of R are called roots and the elements of i?_|_ are called positive 
roots. 

Let GL(V) be the group of invertible linear transformations of V. The Weyl group 
W of the root system R is the subgroup of GL{V) which is generated by the reflections s^ 
(ct G R). The group W is finite and it acts on R. It permutes the possible choices of i?_|_ 
in a simply transitive way. 



3.2. Dunkl operators associated with R. Let i? be a root system in V. Let k:a\-^ 
ka'.R -^ [0,00) be a function which is T^-invariant , i.e., which satisfies ^^,0, = k^ for aU 
w eW and aU a E R. If i? is an irreducible (reduced) root system then the Weyl group 
is transitive on all roots of equal length and there are at most two diff^erent root lengths. 
Thus ka then assumes at most two different values. See the above examples: one root 
length in Ai and two root lengths in C2. The function k is called a multiplicity function. 
The reason for this name is that root systems naturally arise in the structure theory of real 
semisimple Lie algebras, where roots have an interpretation as joint eigenvalues of certain 
operators and the k^i then are (integer) multiplicities of such eigenvalues. 

For ^ E V we will denote by d^ the corresponding directional derivative. The Dunkl 
operators [8], [16] associated with the root system R and the multiplicity function k are 
defined as the operators D^^'"^ : C°°{V) -^ C°°{V) {( G V) given by 

(I?f V)(-) := (5,/)(^) + Y. ^'^K^)^^^?!^- (3-1) 

This definition is easily seen to be independent of the choice of i?_|_. In case of root system 
Ai formula (3.1) reduces for ^ := 1 to formula (2.6). Note that the operator (3.1) consists 
of a term involving a first order derivative and terms involving refiection operators, just as 
we have seen in (2.6). It is an amazing fact, which can be proved in a rather straightforward 

(k) 

way, that the operators Dl commute: 

Let W''' be the algebra generated by the operators D^ . This is a commutative algebra. 

It can be shown that each T^-invariant operator D in W'^' , when restricted in its action 
to the ly-invariant C°° functions on V, coincides with a partial differential operator (so 
its reflection terms vanish when acting on a T^-invariant function). The joint T^-invariant 
eigenfunctions of the T^-invariant operators in W^ are called Bessel functions associated 
with R. In the example Ai things reduce to the one- variable considerations of §2.1 and 
§2.4. More generally, one may study the joint eigenfunctions of the full algebra W''' and 
one may try to do harmonic analysis for these eigenfunctions. A lot of satisfactory results 
have been obtained, see [16] and the references given there. 

3.3. Weight lattice asscociated with R. We still assume a root system R in V. The 
weight lattice P of i? is defined [15] by 

P --{XeV \ 7 ' ; G Z for all a G R}. 

(a, a) 

The subset P_|_ of dominant weights is then given by 

P+ := {A G P I ii^l^ > for aU a G R+}. 



It is easily seen that w(P) = P for w G W, so the Weyl group acts on P. Moreover, it can 
be shown that each Weyl group orbit in P has a one-point intersection with P+: 

VA e P Card(WA n P+) = 1. 

Thus the dominant weights can be used as a set of representatives for the T^-orbits in P. 



P. 
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A:=(4,2) 
Fig. 4. Root system C2 with dominant weights and the set {/j, E P-f | /u -<; A} 

We introduce a partial ordering on P which is induced by the root system: for A, /^ G P 
we say that jj, ~< X ii X — jj, = X^agfl i^a^t for certain nonnegative integers mQ. For root 
system C2 the concepts of this subsection are illustrated in Fig. 4. 

3.4. Trigonometric polynomials associated with R. Let P be the weight lattice of 
a root system R in V. For A G P define the function e on V by 

e^{x) := e*<^'^> {x G V). 

Note that e^ e^ = e^+''. Thus the space 

A := Spanje^ I A G P} 

is an algebra: the algebra of trigonometric functions on V (with respect to R). For a 
function f on V write iwf){x) := f{w-'^x) {w eW, x E V). Then we^ = e""^ {w G W, 
X G P). Put 

mx := J2 e^ (A G P+). 
newx 

Then the functions mx are T^-invariant and they form a basis of the space A^ of W- 
invariant elements in A. 

Let the dual root lattice Q~ be defined by 

Q"-- {XeV \ (A, /i) G Z for all /i G P}. 

This lattice gives rise to a torus 

T := V/i27rQ-'). 

Let a; H^ i be the natural mapping of V oiito T. Then each function / in .4. actually lives 
on T: f{x) = f{x) for a suitable function / on T. 

In the example Ai we have P = "L, P_|_ = {0, 1,2,.. .}, the algebra .4. is spanned by the 
functions x \-^ e*"^ (n G Z) and the subalgebra .4 by the functions 1 and x \-^ 2 cos{nx) 
(n = 1, 2, . . .). The torus T equals R/(27rZ). 



4. Jacobi polynomials associated with R 

4.1. Definition of Jacobi polynomials for R. Let i? be a root system in V and let 
k: R ^ [0, oo) be a ly-invariant multiplicity function as before. Define a weight function 
6k on T by 

Sk{x):= H \2smi{a,x))f^-. (4.1) 

This definition is independent of the choice of i?_|_. Define an inner product on the linear 
space A by 

{f,9)k:= [ f{x)J{^Sk{x)dx {f,geA). (4.2) 

Jt 

Here dx denotes Lebesgue measure on T, normalized such that the volume of T is equal 
to 1. 

The Jacobi polynomial P^ (cf. [11]) of "degree" A G P_|_ and of "order" k is an 
element of A^ of the form 

^f ^ = Yl ^^'/^ ""/- 

such that c\^\ = 1 and 

(Pf \ m^)fc = if /i G P+ and ^ ^ A. (4.3) 

Instead of (4.3) we can equivalently require that P-J^ satisfies the second order differ- 
ential equation 

(a+ 5] k^ coi{\{a,x))d^) P^^\x) = -{\\ + Y.^^^^k^a) P^^\x). (4.4) 

In the example Ai we obtain that Pn \x) = const. C„( cos a;), where C„ is the ultra- 
spherical polynomial of §2.2. The case of the (non-reduced) root system BCi would have 
given us, more generally, the classical one-variable Jacobi polynomials. 

4.2. Three problems and their solutions. As soon as the above definition of Jacobi 
polynomials associated with R is given, three highly nontrivial questions can naturally be 
posed: 

1. It follows immediately from the definition that the orthogonality 

(^f,^i'^)fc = (4.5) 

holds if /^ -< A or A -< ;U. What about (4.5) if A and fj, are not related in the partial 
ordering? 



Prove the existence of a commutative algebra of differential operators with d alge- 
braically independent generators, such that the operators in this algebra have the 
P^ (A G P_|_) as joint eigenfunctions. (Note that the operator in (4.4) can be taken 
as one of the generators.) 

Give an explicit expression for {P^ tP\ )k, or rather for its two factors 



/■pW p('=)\ 

/r>(k) p(fc)\ 
\^0 '^0 /k 



Jt 



and / 6k{x)dx. (4-6) 



In the past few years all these questions have been answered in the positive sense. Let me 
give some indications. 

• If problem 2 can be solved then the answer to 1 follows readily, cf. [11]. Indeed, we 

need sufficiently many differential operators having the P^ as eigenfunctions such 
that the joint eigenvalues, in their dependence on A, separate the points of P_|_. 

• For certain special choices of k the functions P^ , renormalized such that P^ (0) = 1, 
have an interpretation as spherical functions on compact symmetric spaces G/K, cf. 

[14]. (For instance, in case Ai the ultraspherical polynomial Cn can be inter- 

preted as spherical function on the (m — l)-dimensional sphere SO(m)/SO(m — 1).) 
Then problems 1, 2 and the first half of problem 3 can be solved by using the group 
theoretic interpretation. The orthogonality (4.5) for general A, /j, follows by Schur's 
orthogonality relations for matrix elements of irreducible unitary representations of 
G. The first expression in (4.6) was explicitly computed by Vretare [33] in terms of 
Harish- Chandra's c- function related to the spherical functions on the corresponding 
non-compact symmetric space. The algebra of differential operators in problem 2 can 
be obtained by taking the radial parts of the G-invariant differential operators on 
G/K. 

• For the classical root systems question 2 could be answered in a positive way by 
giving explicit expressions for generators of the algebra, see [17] for BC2 and A2, and 
Olshanetsky & Perelomov [26], Sekiguchi [31] and Debiard [6] for the higher rank 
cases. 

• Heckman and Opdam [11] have given positive answers to 2, and hence to 1, by use of 
deep transcendental arguments. This also solved part of Problem 3 (the first expres- 
sion in (4.6)). In 1982 Macdonald [21] had already given conjectures for the explicit 
evaluation of the second expression in (4.6), which could be proved in a number of 
special cases. 

• Problem 3 for general A was solved by Opdam [28] by using so-called shift operators 
[27]. The most simple example, for case Ai, of such operators is the following pair of 
differential recurrence relations for Gegenbauer polynomials: 

^C'^{x) = const. C'^+l{x), 

(1 - x^)-''+^ -^ o (1 - x^)''+^] C^tlix) = const. C'^ix). 
dx I 



iU 



By use of these two formulas we can write J_^{C!^{x))^ (1 — x"^)^ 2 dx as an explicit 

constant times J_i{C^'^i{x))^ {l — x^)'''^^ dx. Opdam's shift operators in general have 
a similar structure of lowering A and raising k, or conversely. The case of root system 
BC2 was already considered in [17], [32]. 

4.3. Dunkl type operators. Some years after Heckman first solved the problems 1 
and 2 of the previous subsection he discovered a dramatical simplification [12] for proving 
these results. For a given root system i? in T^ and a given multiplicity function k he wrote 
down a trigonometric variant of the Dunkl operators (3.1) for ^ E V: 



{D['^f){x):={d^f){x) + l Yl k^W,Ocoti^{c^,x))ifix)-f{s^x)) 

{xeV, /GC°°(y)). (4.7) 

(k) 

Now the operators D^ will no longer commute, in general. However, Heckman showed 

that the operators ^^^vFfl-^^ )"' (^ ^ ^) J = 0)1)2,...), when restricted to the W- 
invariant C°° functions on V, coincide with differental operators which commute with 
each other and form a commutative algebra. This is the algebra looked for in Problem 2 

(k) 

of the previous subsection. The Jacobi polynomials P^ are the joint eigenfunctions of 
the operators in this algebra. By this approach, Heckman also obtained a quick existence 
proof for Opdam's shift operators. 

Next Cherednik [3] made a slight but significant variation in Heckman's Dunkl type 
operators (4.7). He put 

(5fV)(a^):=(<9c/)(^)+ E kc.{c^,0 ^_^\^.A fi^)-fis.x)) 

-\ E ^° (")^)/(^)- 

aG-R+ 

(Here I took the part of the right hand side on the second line from Opdam [13, p. 86]; 
Cherednik is not very specific about this part of his formula.) Cherednik's operators have 
the nice property that they mutually commute, without the need of first restricting to W- 
invariant functions. On the other hand, they do not share the property w D\ w~^ = D^J 
of Heckman's operators. Anyhow, by means of Cherednik's operators one can draw the 
same conclusions as by Heckman's operators, and in a similar way. Moreover, a structure 
of graded Hecke algebra can be associated with Cherednik's operators. 



ii 



5. Macdonald polynomials associated with R 

5.1. Definition of Macdonald polynomials. Let < q < 1. We keep the assumptions 
of §4.1 except that we replace the weight function 6k in (4.1) by 



^k,q{x) := Y[ 
cteR^ 



(e*<"'^^;g)c ^ 



/CXD 



(^gfccg*(a,a:>- gr)^ 



Then the Macdonald polynomials P^ ' were defined by Macdonald [23], [24] just as the 
Jacobi polynomials P^ , but with the inner product in (4.2) replaced by 



f,g)k,q-= / f{x)g{x)Sk,q{x)dx {f,geA). 
Jt 



In the case of root system Ai the Macdonald polynomials coincide with the g-ultraspherical 
polynomials x h^ C^'^(cosa;). For any root system R, in the limit for q ^ 1, the Macdonald 

(kg) ... (k) 

polynomial P^ ' tends to the corresponding Jacobi polynomial P^ . 

(kg). 

Macdonald gives some explicit g-difference operators of which the P^ ' are eigenfunc- 
tions. Although these operators, except for root system An (where they were independently 
found by Ruijsenaars [30]) do not yet give a full commutative algebra of operators having 
the P^ ''^ as joint eigenfunctions, the additional parameter q gives enough freedom such 
that already the eigenvalue of one such operator separates the elements of P_|_ for generic 
q, by which a positive answer to question 1 in §4.2 can be given for the case of Macdonald 
polynomials. Taking limits for q ^ 1 then yields the same positive answer for the case of Ja- 
cobi polynomials. This is an alternative to Heckman's approach via Problem 2. Macdonald 
also gives conjectured explicit expressions for the squared norms (P^ '"^ , P^^ '"^ )k,q- 

5.2. Askey- Wilson polynomials for root system BCn- The author [18] introduced 
for the non-reduced root system BCn a class of polynomials having two more parameters 
than Macdonald's class for BCn- This extended class reduces for n = 1 to the Askey- Wilson 
polynomials [2]. In [18] only one explicit g-difference operator was given having the BCn- 
polynomials as eigenfunctions, but this was sufficient for establishing orthogonality. Later, 
van Diejen [7] gave explicit expressions for the generators af a full commutative algebra of 
operators having the BCn polynomials as joint eigenfunctions. 

5.3. Cherednik's approach to Macdonald polynomials. Cherednik [4], [5] suc- 
ceeded to give positive answers to questions 2 and 3 in §4.2. In the context of certain 
representations of affine Hecke algebras he could realize a commutative algebra of opera- 
tors which have the Macdonald polynomials as joint eigenfunctions. In the same context 
he could realize g-analogues of Opdam's shift operators and next, by the same technique 
as in Opdam, prove Macdonald's conjectures in the Q-case. 

It is beyond the scope of this short survey to explain Cherednik's approach in any 
detail. In May 1994 I. G. Macdonald delivered some very helpful lectures in Leiden in order 



to explain Cherednik's approach. Let me here only give a few indications. Just as a Hecke 
algebra is a deformation of the group algebra of a Weyl group, an affine Hecke algebra (cf. 
[20]) deforms the group algebra of an affine Weyl group. If R is an irreducible root system 
in V with Weyl group W then the (extended) afRne Weyl group is the seniidirect product 
W :=W \K P~, where the dual weight lattice P^is defined as P' := {X e V \ (A, a) € Z}, an 
abelian group under addition. Then W acts as a group of motions on V, with P" acting as 
a group of translations. The group W also acts on A, with W acting as before and with 
the action of P~' still depending on a parameter q. 

The afBne Hecke algebra H can be defined in terms of generators and relations which 
still depend on the values of a T^-invariant function a h^ t^ on R. Corresponding to a 
choice of i?_|_ we can define Pj^ . The embedding of Pj^ in H then generates a commutative 
subalgebra y oi H . 

For given q we can use the action of T^ on ^ in order to define an action of H on A, 
by specifying the action for a set of generators of H (Demazure operators). This action 
depends on q and the t^. Put to, = g~'^«/2^ where k: R ^ [0, oo) is a multiplicity function. 

(kg) 

Then Cherednik proves that the Macdonald polynomials P^ ' are the joint eigenfunctions 
of the T^-invariant elements in the commutative algebra y. This answers question 2 in 
§4.2. 

At the moment it is still an open problem to extend Cherednik's approach to the BCn 
polynomials. As van Diejen [7] already answered question 2 in §4.2 in a positive way for 
this case, it would be nice to complement van Diejen's constructive approach with the deep 
conceptual approach via affine Hecke algebras. 
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